Nonequilibrium characteristics in all-superconducting tunnel structures 



M. A. Laakso, 1 P. Virtanen, 1 F. Giazotto, 2 and T. T. Heikkila^B 

1 Low Temperature Laboratory, Helsinki University of Technology, P.O. Box 2200 FIN-02015 TKK, Finland 
2 NEST-CNR INFM & Scuola Normale Superiore, 1-56126 Pisa, Italy 
(Dated: February 6, 2008) 

We study the nonequilibrium characteristics of superconducting tunnel structures in the case when 
one of the superconductors is a small island confined between large superconductors. The state of 
this island can be probed for example via the supercurrent flowing through it. We study both the 
far-from-equilibrium limit when the rate of injection for the electrons into the island exceeds the 
energy relaxation inside it, and the quasiequilibrium limit when the electrons equilibrate between 
themselves. We also address the crossover between these limits employing the collision integral 
derived for the superconducting case. The clearest signatures of the nonequilibrium limit are the 
anomalous heating effects seen as a supercurrent suppression at low voltages, and the hysteresis 
at voltages close to the gap edge 2A/e, resulting from the peculiar form of the nonequilibrium 
distribution function. 



I. INTRODUCTION 

New device concepts based on nonequilibrium ef- 
fects in superconducting mesoscopic tunnel structures 
have been proposed in the last few years. These in- 
clude Josephson transistors, electron refrigerators and 
thermometers J« In Josephson transistors the supercur- 
rent flowing through a superconductor-normal metal- 
superconductor (SNS) weak link can be suppressed or 
even reversed in a 7r-transition2id£ by driving the normal 
metal part out of equilibrium through injection of charge 
carriers from additional terminals. When the additional 
terminals are superconductors connected by tunnel junc- 
tions, the supercurrent can also be enhanced^ This tran- 
sistorlike operation with large current and power gain 
has also been experimentally demonstrated.— Also an 
all-superconducting SISIS transistor in the quasiequilib- 
rium regime has been theoretically addressed^ In the 
quasiequilibrium limit the electron-phonon interaction is 
nearly absent and the sample can be considered as de- 
tached from the phonon bath. The high frequency of 
electron-electron collisions still serves as a method of re- 
laxation and the electrons assume a Fermi distribution 
but with a temperature that in general differs from the 
temperature of the phonon bath. Here we study a simi- 
lar SISIS structure with arbitrary strength of the inelas- 
tic scattering seeking ways to characterize the degree of 
nonequilibrium on the system. The paper is organized as 
follows: The model of the SISIS structure is presented in 
Sec. 2. All the relevant equations and calculated results 
are presented in Sees. 3 and 4, respectively. We finish 
with a summary and a discussion in Sec. 5, where we 
also address briefly the feasibility of this structure as a 
transistor. 

II. MODEL 

The superconducting structure under study is 
schematically depicted in Fig. [TJ We characterize the 
mean free path that the electron travels before scatter- 



ing by scattering lengths l e i for elastic scattering and 
l e ~ P h and l e -e for inelastic electron-phonon and electron- 
electron scattering, respectively. In mesoscopic systems 
typical orders of magnitude are l e i « 10 . . . 100 nm and 
l e - e ~ 1 ■ ■ ■ 20/xm. The electron-phonon scattering length 
depends strongly on temperature. For a typical copper 
wire l e -ph ~ 21/xm at 1 K but at 100 mK we already have 
le-ph ~ 670 jiia— In other metals these length scales are 
of the same order of magnitude. The superconducting is- 
land in the middle is assumed to have small dimensions so 
that L <C l e _ ei l e _ p f l , leading to weak energy relaxation 
via inelastic scattering. As we shall see in the following, 
in this case it is possible to drive the electron energy dis- 
tribution of the island out of equilibrium by quasiparti- 
cle injection from the superconducting leads. The degree 
of nonequilibrium on the island can then be probed for 
instance by measuring the supercurrent driven through 
the island via an additional SISIS line. The leads are 
assumed to remain in thermal equilibrium due to their 
large dimensions. We further assume that the resistances 
of the tunnel contacts are large compared to the normal 
state resistance of the superconducting island. This al- 
lows us to use the tunnel Hamiltonian approach, in which 
each region has spatially constant, separate energy distri- 
butions, independent of the directions of the momenta. 



III. FORMALISM 

A. Green's functions in SISIS structure 

We use the quasiclassical Keldysh Green function for- 
malism together with the tunnel Hamiltonian model in 
describing our system. It has previously been successfully 
applied to hybrid structures with normal metal and su- 
perconducting islands The quasiclassical Green func- 
tions in Nambu space can be written in a matrix form 
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FIG. 1: Scheme of the SISIS structure studied in this work. 
The superconducting island (2) in the middle is connected 
with tunnel contacts to four large superconducting leads 
(1,3,4 and 5). The control line is biased with voltage V, which 
controls the energy distribution on the island. A supercurrent 
Is is driven across the island from lead 4 to lead 5, and its 
magnitude depends on the distribution function on the island. 
Each SIS junction is a tunnel contact of resistance Ri. 



Here A is the superconducting order parameter, 7 = 

J2j Vjdf an d S = J2j Vjff^- The index j runs over 
the other parts of the structure that are connected with 
tunnel contacts to the region in question. The character- 
istic tunneling rate 77 between superconductors is defined 
as rjj — (4i/e 2 flRj)~ 1 , where v is the normal-state den- 
sity of states at the Fermi level and fl is the volume of 
the island. In the tunneling limit jj <C A and we may 
neglect the exact forms of 7 and S and instead use some 
constant 7 and S in the numerical simulations. Below, we 
choose 7 = 1CP 4 and 6 = 10~ 5 . This value of 7 has been 
experimentally verified in Ref. 2. The Keldysh Green 
function for the system can be written with the standard 
parametrization as 



_ 9 K = 9 R (fL + hfr) - (h + t 3 /t).9 A , (3) 
where T3 is the third Pauli spin matrix. Here we have 
also used the odd- and even-in- E parts of the distribution 
function: 



f L (E)=-f(E) + f(-E), 
f T (E) =1 - f(-E) - f(E). 



as 



9 f 
~P 9 



(1) 



where g is cither retarded (advanced), g ^ ' , or Keldysh, 
g K , Green's function. In the tunnel Hamiltonian model 
Green's functions are isotropic with respect to directions 
of the momenta. In this case the retarded (advanced) 
functions satisfy the steady-state Usadel equations*^ 
with solutions 
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The full distribution function can be recovered with 
2f(E) = 1 - f L (E) - f T {E). We also define 
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The functions /' + ' and are even in E and f(~^ is 
odd. The density of states is given by g(~'. The odd and 
even parts of the noncquilibrium distribution function 
can now be found from the kinetic equations presented 
in Ref. 9. The resulting equations are 



-AJ lV2 e 2 n 2 =g[ jGj {g[J +fll (hi + hi - hi - hi) + g[ (hi - hi - hi + hi)} 

+ 9i\e G 3 { 03^+^3 (? L2 + hi - fL3 - hi) + 9 { 3 ~E-^ (hi - /T2 - fh3 + hs)} , 

(4) 

(8|A 2 |/ T2 ^ +) - 4J 2 ) u 2 e 2 n 2 = 9 Qg 1 (hi + hi - hi - hi) + sB-^ (~hi + hi + hi ~ hi)} 

+ 9 { 2 ~E G 3 { 03^+^3 (hi + hi ~ fh3 ~ hi) + g^E-^ (~hl + hi + fh3 - fT3)} , 

(5) 

I 

where Gi = 1 /Ri are the conductances of the tunnel and Ji are the collision integrals for the energy relaxation, 
contacts, ^ are the chemical potentials of the regions i 
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B. Order parameter and currents 

The pair potential in the central island must be solved 
self-consistently from the equation 

|A 2 | = -/ dEf L2 f { E \ (6) 

1 J-Ec 

where Eq is the BCS cutoff energy and A is the electron- 
electron interaction parameter. When a SIS-junction is 
not biased with an external voltage, the supercurrent 
flowing across the junction is given by 

x sin( X4 - X2) + (/T2 - hi) (g^oV + A +) A +) 

X C0S(X4 - XiYj } ■ (J) 

The first part of the equation multiplying the sine term 
represents the usual dc-Josephson relation where X2A 
are the macroscopic phases of the respective supercon- 
ductors. The term j\ ' is finite only between 

A2 < E < A4 whereas the term /j ^ is finite when 
A4 < E < A2. The second part in Eq. ([7]) usually van- 
ishes because h = in quasiequilibrium. If a finite 
charge imbalance develops on the island, Jti deviates 
from zero, and the second part contributes to the cur- 
rent as well. 

If a voltage is applied across the junction the phase 
difference begins to evolve in time and the supercurrent 
averages to zero. In this case only the tunneling of quasi- 
particles contributes to the current, so that we have 

1 r°° r 

I 1 ^ 2 =-^j M ^ {ffU+^U UL2 + fT2 

-hi - hi) + 9[:LA: ] e + /t2 

+/li-/ti)}. (8) 

The quasiparticles tunneling through the junction also 
carry heat. The energy current is 

I f°° ( 

^ = ~ I^Rl ^ X^+A^ (/« + h2 

-hi - hi) + 9 ( Ce-^e (h2 - hi 

-hi + hi)}, (9) 

which is used in determining the electron temperature in 
quasiequilibrium. 



can be included with e — e collision integrals. We may 
also study cases where the dimensions of the island are 
no longer significantly smaller than the electron-electron 
scattering length, i.e., L < l e -e <C l e -ph- The collision 
integral for a screened Coulomb interaction in a diffu- 
sive wire is known^ and has been used in the analysis 
of a SINIS system.— It is strictly valid only for a normal 
metal island, however. To get a qualitative picture of the 
changes due to superconductivity in energy relaxation, 
we apply instead a collision integral where the structure 
of the Nambu space has been taken into account. In 
the clean limit the potential of a distant electron is com- 
pletely screened by all other electrons in the superconduc- 
tor and the electron-electron interaction can be approx- 
imated by a point interaction. In this case the potential 
may be modelled with a delta function V(r) = U2X ee S(r) 
and the collision integral is^ 

(E 3 ) =njj dE ld E 2 { - /£>/£>) 

-h (1 - f El ) (1 - h 2 ) (1 - h 3 ))}, (10) 

where k = 4Al e ir/pFVF, vp and pf are the Fermi veloc- 
ity and momentum, respectively, and energies satisfy the 
conservation law E = E\ + E2 + E3 . The second collision 
integral, J<f e \ vanishes in a left-right symmetric struc- 
ture. We note that because the terms o^T and fiT in 
the kernel of the integral assume very small values when 
I-E3I < A, the collision integral has a very small effect 
on excitations inside the gap. On the other hand, energy 
relaxation is strongest for excitations at | £^3 1 = A due to 
sharp peaks at the edge of the gap in these same terms. 



IV. RESULTS 

A. Full nonequilibrium 

We begin by presenting the calculated distribution 
function along with the order parameter and electric 
currents for the simplest, namely left-right symmetric 
case, where the tunnel junction resistances are the same 
and reservoirs 1 and 3 are similar superconductors, i.e., 
Ri = R 3 = R and |Ai| = |A 3 | = |A L |. When the struc- 
ture is biased with a voltage V, the conservation of elec- 
tric current forces the chemical potentials of reservoirs 1 
and 3 to /ii = eV/2 and fi^ = —eV/2, respectively. 



C. Energy relaxation 

In practice the inelastic scattering is never completely 
absent. At low temperatures the most relevant relax- 
ation mechanism is electron-electron scattering, which 



1. Distribution function 

The solution of the kinetic equations ((3]) and (JSJ) in 

the absence of energy relaxation {J\ — J 2 — 0) may be 
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FIG. 2: (Color online) Nonequilibrium distribution function 
for the superconducting island at T — 0.1 Tc- The cooling 
effect reducing the number of excited quasiparticles as the 
voltage is increased is evident. Here and below we denote 
Ao = Al(T — 0) and Tc is the critical temperature of the 
leads. 

written in terms of the full distribution functions as 

, 9E+fJi + 9 { E -y.h , 
h - — — — — ■ i 11 ) 

This form is remarkably simple due to the symmetry of 
the problem, and it can also be derived by considering the 
conservation of electric current^ The distribution func- 
tion is plotted in Fig.[2]for various bias voltages at a bath 
temperature of 0.1 Tc- The critical temperature of the 
superconductor is Tc — (1.764feg) _1 Ao- With a small 
voltage bias fewer of the states below the Fermi level are 
occupied whereas the occupation is increased above the 
Fermi level. This increase in excited quasiparticles can 
be interpreted as a heating of the island. This anoma- 
lous heating effect stems from the assumption of a finite 
7 in Eq. i.e., from the presence of quasiparticle states 
within the gap. In the absence of these states, no anoma- 
lous heating is observed. Once the voltage is increased 
above eV = Al, the number of excited quasiparticles on 
the island begins to decrease due to extraction to states 
right above the energy gap in the superconducting reser- 
voirs. This cooling effect is discussed in Refs. 1 and 2. In 
Fig. [3] the distribution function is plotted at higher bath 
temperatures for a bias voltage eV/A — 3. At these 
temperatures the reservoirs have more excited quasipar- 
ticles above and below the gap, and the small notches at 
\E\ = eV/2 + Al are a result of their injection. 

2. Order parameter 

In order to measure the degree of nonequilibrium on 
the island we must look for nonequilibrium induced ef- 



FIG. 3: (Color online) Nonequilibrium distribution function 
for the superconducting island at various bath temperatures 
for a voltage eV/ Aq — 3. 

fects in some measurable quantities, e.g., supercurrent 
through the island. First we calculate the magnitude 
of the order parameter with the self-consistency equa- 
tion ^ . In general this must be solved numerically. The 
magnitude of the order parameter of the island as a func- 
tion of voltage at various bath temperatures is shown in 
Fig. Ha). At T — 0.1 T c the odd-in-£: part of the distri- 
bution is effectively unchanged outside the gap giving the 
same result as for equilibrium. However, once eV > 2Al 
the peculiar shape of the distribution makes it possible 
to have a lower value solution for the order parameter as 
well, giving rise to a hysteretic behavior with three solu- 
tions. Once voltage reaches eV = 2(A2 + Al) only the 
smallest solution, namely A2 = 0, is possible. This is due 
to the fact that the order parameter can never exceed its 
zero-temperature value. The multivalued behavior of the 
order parameter can be interpreted as different minima 
and maxima in the free energy^i^ In this case the largest 
and smallest values represent minima and the middle 
value represents a maximum. If we increase the volt- 
age from zero, the system stays in the free-energy min- 
imum corresponding to a superconducting state. Once 
we enter the hysteretic region, thermal fluctuations may 
cause the system to jump to normal state, which is the 
other free-energy minimum. In the absence of fluctua- 
tions, the system finally jumps to the normal state at 
eV = 2(A 2 + Al). If we now proceed by decreasing the 
voltage, the jump to the superconducting state may again 
occur somewhere in the hysteretic region. Once the volt- 
age is decreased enough, only the superconducting state 
is possible. 

At higher bath temperatures the order parameter is 
initially in its equilibrium value, but increases along 
the voltage as the island cools. In Fig. |4jb) the or- 
der parameter is shown at T = 0.7 Tc but for different 
zero-temperature ratios A 2 /A L = T C2 /T C - For ratios 
A2/ Al < 0.7, the island is initially in the normal state 
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FIG. 5: (Color online) Supercurrent through the island in 
full nonequilibrium as a function of bias voltage at various 
bath temperatures with a ratio A2/A.L = 0.3. The arrows 
indicate the direction the curve is traced when the bias voltage 
is varied. Thermal fluctuations may cause the discontinuous 
jump to occur somewhere in between the two extremes shown 
in the figure. The system is assumed symmetric, i.e., Ri = 
R3 = R4. = Rs and A 4 = A 5 = A L . 



FIG. 4: (Color online) Order parameter as a function of bias 
voltage at various bath temperatures (a) and ratios A2/A1, 
(b). 



because the bath temperature is above its critical tem- 
perature. Upon increasing the voltage, the island turns 
superconducting once the electron distribution has fea- 
tures sharp enough to support an energy gap. 



3. Electric currents 

Now we examine the effect that the magnitude of 
the order parameter has on the electric current driven 
through the island. In light of the results in the pre- 
vious subsection the measurements should be made at 
relatively high temperature in order to fully bring out 
the variation in the energy gap. By choosing a setup 
with a lower {±.2/ ratio enables us to use a lower abso- 
lute temperature and thereby also minimize the electron- 
phonon relaxation, because the power injected into the 
phonons depends on temperature as T 5 J^ The supercur- 
rent through the island is calculated with Eq. and it 
is presented in Fig. [5] for various temperatures assuming 
a ratio A2/A^ = 0.3, which corresponds roughly to the 
Ti/Al combination. When the bath temperature is lower 
than the critical temperature of the island, the initial 



heating effect with low bias voltages is evident. In the 
cooling regime the bath temperature has a negligible ef- 
fect on the magnitude of the supercurrent. The hysteresis 
of the order parameter carries over to the supercurrent 
but no 7r-state is observed. It is illustrative to compare 
these to the corresponding results in quasiequilibrium, 
where the high frequency of electron-electron collisions 
force the quasiparticlcs on the island to assume a Fermi 
distribution. The electron temperature in quasiequilib- 
rium can be obtained by demanding that the energy cur- 
rent in Eq. © to the island vanishes (we also neglect the 
contribution of electron-phonon interaction to the energy 
current). The supercurrent in quasiequilibrium is shown 
in Fig. [5] In quasiequilibrium the heating effect is ab- 
sent and the island cools even with low voltages resulting 
in an increase of the supercurrent. Superconductivity is 
lost once the voltage exceeds eV — 2(A 2 -I- A^), just as 
in full nonequilibrium. The falling edge here is not hys- 
teretic, however. A further means to probe the degree 
of nonequilibrium is to voltage bias the second SISIS-line 
as well and measure the energy gap from the I-V curve. 
The quasiparticle current flowing through the probe junc- 
tion in this case may be calculated with Eq. ([H]). The 
resulting I-V curve does not differ from its equilibrium 
shape, in which the current has a discontinuos jump at 
eV = 2(A 2 + A L )^ The value of A 2 and its hystere- 
sis changes the voltage at which the jump is observed, 
however. 
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FIG. 6: (Color online) Supercurrent through the island in 
quasiequilibrium as a function of bias voltage with parameters 
identical to the full nonequilibrium case presented in Fig. [5] 



current at voltages over eV = 2. This corresponds to a 
transition above which A4 > A2. By choosing a setup 
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FIG. 7: (Color online) Distribution function for the supercon- 
ducting island at eV/Ao = 1.5 (left) and eV/Ao = 3 (right) 
for various IC co u with T = 0.7 Tc- K, C oll = 00 corresponds to 
quasiequilibrium . 



B. Nonequilibrium with energy relaxation 

When the energy relaxation due to inelastic electron- 
electron scattering is taken into account, we are no longer 
able to obtain an explicit expression for fa. In the left- 
right symmetric case we must instead solve the resulting 
integral equation 

_ 9 ( E^h+9 ( E ^h + (e 2 ^ 2 iWg<:i) j[ ee) [f 2 ] 

The relaxation strength can be adjusted by varying the 
parameter JC co u — fte 2 i/2^2-Ri- The distribution func- 
tion calculated for various values of K, co u is shown in 
Fig. [7l The energy distribution gradually relaxes to- 
wards a Fermi distribution upon increasing the strength 
of the relaxation. The influence of inelastic scattering 
to the supercurrent is shown in Fig. [8] for a structure 
consisting entirely of one type of a superconductor. The 
enhancement of superconductivity is suppressed as the 
electron-electron collisions drive the electron tempera- 
ture of the central island towards quasiequilibrium. With 
the strongest relaxation the cooling effect is completely 
lost and the supercurrent drops smoothly to zero as the 
voltage is increased. With the two largest strengths of re- 
laxation the hysteresis is lost as well. At larger voltages 
the supercurrent is a non-monotonous function of K, co u , 
as the supercurrent in quasiequilibrium (JC co u = 00) is 
significantly larger than with moderate relaxation. This 
can also be seen in the left part of Fig. \7\ where the 
distribution function in quasiequilibrium is sharper com- 
pared to the distribution with IC co u = 10. This sharpness 
leads directly to a larger supercurrent. The curves with 
ICcoii = 1 and IC C oii = 2 show a small jump in the super- 



where leads 4 and 5 have a different energy gap from the 
rest of the system, this peak could be seen at different 
voltages. 
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FIG. 8: (Color online) Supercurrent through the island as a 
function of voltage for various K, co ii with T — 0.7 Tc- For the 
hysteretic curves the arrows indicate the direction the curve 
is traced when the bias voltage is varied. 



C. Asymmetric structure 

Let us now examine an asymmetric situation, where 
Ri ^ i?3 or Ai 7^ A3. By solving the kinetic equations 
(H]) and (J5J) without relaxation we obtain quite lengthy 
expressions for the odd and even parts of the distribution 
function 
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FIG. 9: (Color online) Supercurrent through the island as a 
function of voltage for different degrees of asymmetry in the 
SISIS control line. The arrows indicate the direction the curve 
is traced when the bias voltage is varied. The inset shows the 
distribution function on the island for a ratio R1/R3 = 2 
at eV/Ao = 1.5. The distribution function exhibits small 
asymmetry due to finite fx as can be seen by the additional 
notch at negative energies. 



The distribution functions depend on the volume, energy 
gap and normal state density of states at the Fermi level, 
but these can be included in dimensionless constants of 
the type G/\A\vfle 2 . In the asymmetric case the poten- 
tials /ii and /Z3 must be chosen such that the electrical 
current is conserved. This implies the vanishing of the 
total net current into the island, i.e., I 1 ^ 2 = j 2 ^ 3 cal- 
culated with Eq. ©. 

The supercurrent for xi ~ Xi — f/2 and different 
ratios Ri/R 3 is shown in Fig. [5J In an asymmetric 



FIG. 10: (Color online) Potential difference between quasipar- 
ticles and the condensate for different degrees of asymmetry 
in the SISIS control line. 



structure the magnitude of the order parameter seems 
to be close to its value in the symmetric case with a 
voltage of eV — 2max(|/xi|, I/K3Q. This is reasonable 
because the distribution function in the region \E\ > 
max(Ai + I , A 3 + |/x 3 |) is similar to the distribution 
in the symmetric structure as shown in the inset. Su- 
perconductivity is lost once > A 2 + Ax or |/i 3 | > 
A 2 + A3. With high asymmetry ratios the potentials 
differ very much from ±eV"/2 and superconductivity is 
lost at a lower bias voltage compared to the symmetric 
structure. Also the hysteretic region is evident. 

Because the charge imbalance function, fx, is finite, 
also the latter part of Eq. (JT]) may contribute, depend- 
ing on the phase difference between the superconductors. 
Its magnitude can be investigated by setting the phase 
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difference to zero. In this case the electric current is 
significantly smaller, of the order of 10~ 3 eR^I/Ao, and 
mostly due to quasiparticle current induced by the charge 
imbalance. The charge imbalance leads to a difference in 
the chemical potentials between quasiparticles and the 
condensate. The potential difference is given by^ 

e* = - PJ^fagV. (14) 

This quantity is shown in Fig. [10] for the superconduct- 
ing regime. The quasiparticle current depends linearily 
on this potential difference. The equation for the super- 
current also seems to imply a cos(Ax) dependence in the 
supercurrent. This deviation from the dc Josephson re- 
lation is negligible however, because the integral over the 
supercurrent term fi^ is three orders of magnitude 
smaller than over the quasiparticle current term ^ g\ \ 

V. DISCUSSION 

According to our results there are several mea- 
surable features present in a nonequilibrium, all- 
superconducting, tunnel structure. The initial electron 
heating is seen as a strong suppression in superconduc- 
tivity of the central island when the tunnel structure is bi- 
ased with a low voltage. This is observable when the bath 
temperature is slightly below the critical temperature of 
the central island but well below the critical tempera- 
ture of the superconducting leads. The nonequilibrium 
cooling effect together with the destruction of supercon- 



ductivity at eV = 2(A 2 + A^) should be observable with 
a wide range of configurations. The accompanying hys- 
teresis with low or nonexistent relaxation can be seen in 
the supercurrent as well. The magnitude of the energy 
gap could be directly measured with a quasiparticle cur- 
rent probe, where the jump in the current happens at a 
probe voltage of eV = 2(A 2 + A^). 

Due to hysteresis the application of this structure as 
a transistor is unfeasible in states far from equilibrium. 
With moderate to strong relaxation the hysteresis is ab- 
sent and does not hamper transistor-like operation. The 
sharp current-voltage characteristics giving rise to high 
differential current gain are lost with the strongest relax- 
ation, however. If actual power gain were to be achieved, 
the Josephson junctions have to be operated in the dis- 
sipative regime and coupling to the environment should 
be taken into account in the calculations. 

Small asymmetries of the order of ten percent in the 
system do not give rise to qualitatively different be- 
haviour. Asymmetries larger than that begin to develop 
charge imbalance in the central island, leading to differ- 
ent chemical potentials for the superconducting conden- 
sate and quasiparticle excitations. This potential differ- 
ence can be observed in the quasiparticle current flowing 
to the island from both reservoirs 4 and 5, when the phase 
difference accross the Josephson junctions vanishes. 
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